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Abstract 

We study the propagation of gravitational waves carrying arbitrary 
information through isotropic cosmologies. The waves are modelled 
as small perturbations of the background Robertson- Walker geometry. 
The perfect fluid matter distribution of the isotropic background is, 
in general, modified by small anisotropic stresses. For pure gravity 
waves, in which the perturbed Weyl tensor is radiative (i.e. type N 
in the Petrov classification), we construct explicit examples for which 
the presence of the anisotropic stress is shown to be essential and 
the histories of the wave-fronts in the background Robertson- Walker 
geometry are shear-free null hypersurfaces. The examples derived 
in this case are analogous to the Bateman waves of electromagnetic 
theory. 
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1 Introduction 



This paper is primarily concerned with investigating the influence of grav- 
itational waves on the matter content of the universe. The gravitational 
waves are modelled as small perturbations of Friedmann-Lemaitre cosmolog- 
ical models with a spatially homogeneous and isotropic Robertson-Walker 
(RW) geometry We use the gauge- invariant and covariant perturbation 
approach of Ellis and Bruni ]IJ. This technique has previously been used 
to study aspects of gravitational wave propagation in isotropic cosmologies 
which differ from those considered here (see, for example, [@], 0). 

We look for gravitational waves which carry arbitrary information and 
this motivates us to require the Ellis-Bruni gauge-invariant variables to have 
an arbitrary dependence on a function. Thus by "arbitrary information" in 
the perturbations we mean that they depend upon an arbitrary function. 
This in turn means that the profile of the waves described by the perturba- 
tions is unspecified. The idea of introducing arbitrary functions into solutions 
of Einstein's equations describing gravitational waves goes back to pioneering 
work by Trautman M. This point of view was initiated in the present con- 
text by Hogan and Ellis H, but whereas in that work the perturbed matter 
distribution was a perfect fluid we here allow it to be completely general. 
The Ellis-Bruni approach does not involve working directly with Einstein's 
field equations but instead the basic equations used are the Ricci identities, 
Bianchi identities and the matter equations of motion and energy conser- 
vation equation, with Einstein's field equations incorporated in them. We 
work in the linear approximation in terms of the Ellis-Bruni variables and 
demonstrate that consistency of our assumptions with the basic equations 
necessarily leads to all perturbation variables vanishing except the perturbed 
shear of the matter world-lines and the anisotropic stress perturbation of the 
matter distribution. The consistency of the equations satisfied by these sur- 
viving variables is established. We then specialise our study to pure gravity 
wave perturbations for which the perturbed Weyl tensor is radiative (type 
N in the Petrov classification). We show that for such waves the histories 
of their wave-fronts in the background RW space-time can be shear-free 
null hypersurfaces. The various open and closed RW geometries admit, in 
a natural way, families of shear-free null hypersurfaces. We derive explicit 
solutions of the perturbation equations describing gravitational waves propa- 
gating through these isotropic universes for which the histories of their wave 
fronts are these naturally occuring families of shear-free null hypersurfaces. 
These solutions are analogous to the Bateman waves in electromagnetic 
theory. 

The paper is organised as follows: in section 2 we list the basic exact 
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equations derived from the Ricci identities, Bianchi identities and the equa- 
tions of motion and of energy conservation of the matter distribution. In 
section 3 we introduce the Ellis-Bruni gauge-invariant variables and show 
how for us they have an arbitrary dependence on a function. The consistency 
of these assumptions with the basic equations in the linear approximation is 
then systematically studied, and the mathematical consistency of the surviv- 
ing equations is established. This is followed in section 4 by specialisation 
to pure gravity wave perturbations and the demonstration that their wave- 
front histories can be shear-free null hypersurfaces. In section 5 some explicit 
families of pure gravity wave perturbations are derived. The paper ends with 
a discussion in which we comment on the solutions obtained in section 5. 



2 The Basic Exact Equations 

To make the paper as self-contained as possible the basic equations re- 
quired for our study are given in this section. We use the notation and sign 
conventions of |7| and all of the equations given here can be found in 0] 
with the exception that the Bianchi identities given in J7j apply to a perfect 
fluid matter distribution whereas we require the extension of these to a gen- 



eral matter distribution with energy-momentum-stress tensor given by ( 2.2 
below. This covariant approach to cosmology began in a systematic way 
with the work of Schucking, Ehlers and Sachs [see [[?], [Q for example] while 
Hawking || gave the first description of cosmological perturbations in this 
context. We are concerned here with a four dimensional space-time manifold 
with metric tensor components g ab , in a local coordinate system {x a }, and 
a preferred congruence of world-lines tangent to the unit time-like vector 
field with components u a (with u a u a = —1). With respect to this 4-velocity 
field the Weyl tensor, with components C abcd , is decomposed into an "electric 
part" and a "magnetic part" given respectively by 

E ab = C apbq vP u" , H ab = *C apbq vP u q . (2.1) 

Here *C apbq = ^rj ap rs C Tsbq is the dual of the Weyl tensor (the left and right 
duals being equal), rj abcd = sf=g e abcd with g = det (g ab ) and e abcd is the Levi- 
Civita permutation symbol. The expression for the Weyl tensor in terms of 
E ab and H ab is given in . The symmetric energy-momentum-stress tensor 
T ab can be decomposed with respect to the 4-velocity field u a as 

T ab = fiu a u b + ph ab + q a u b + q b u a + rr ab , (2.2) 

where 

h ab = g ab + u a u b , (2.3) 
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is the projection tensor and 



q a u a = 0, 7r ab u b = 0, 7r a a = 0, (2.4) 

with ir ab = n ba . Here \x is interpreted as the matter energy density mea- 
sured by the observer with 4-velocity u a , p is the isotropic pressure, q a is 
the energy flow (such as heat flow) measured by this observer and n ab is 
the anisotropic stress (due, for example, to viscosity). We shall indicate co- 
variant differentiation with a semi-colon, partial differentiation by a comma, 
covariant differentiation in the direction of u a by a dot and a definition by a 
colon followed by an equality sign. Thus the 4-acceleration of the time-like 
congruence is 

u a :=u a ;b u b , (2.5) 
and u a -b can be decomposed into 

u a; b = uj ab + a ab + ^0 h ab -u a u b , (2.6) 

with 

UJab ■= U[ a - b ] + U[ a U b ] , (2.7) 

the vorticity tensor of the congruence tangent to u a . The square brackets 
denote skew-symmetrisation as usual. Also 

<?ab ■= U(o;6) + «(o «6) - \@ K b , (2.8) 

is the shear tensor of the congruence. The round brackets denote symmetri- 
sation and 

9 ■= u a , a , (2.9) 

is the expansion (or contraction) of the congruence. 

The key equations we shall need are obtained by projections in the direc- 
tion u a and orthogonal to u a (using the projection tensor h ab ) of the Ricci 
identities 

Ua;dc - Ua-cd = RabcdU b , ( 2 -10) 

where R a bcd is the Riemann curvature tensor, the equations of motion and 
the energy conservation equation contained in 

T ab ;b = 0, (2.11) 

and the Bianchi identities written conveniently in the form 



Qabcd j^c[a;b] _ c[a 



■ 4 = ^ ---g c[a K^ . (2.12) 



4 



Here R ca := R cba b are the components of the Ricci tensor and R := R c c is the 
Ricci scalar. Einstein's field equations, after absorbing the coupling constant 
into the energy-momentum-stress tensor, take the form 

Rab — ^ 9ab R = T ab . (2-13) 
The Ricci identities yield Raychaudhuri's equation, 

6+ 1 -6 2 -u? a + 2(o- 2 -uj 2 ) + 1 -^ + 3p) = } (2.14) 
(here a 2 := \a ab a ab , uj 2 := ^oj ab uj ab ) the vorticity propagation equation, 

h a b u b + ^6u a = a\ u b + ^ V abcd u b ii c . d , (2.15) 

(here tu a := \ rf bcd u b u cd is the vorticity vector) the shear propagation equa- 
tion, 

2 

K K (<jfg - u (f;g)) ~u a u b + uj a uj b + a a} a f b + - 9 a ab 

+h ab (-^co 2 - -a 2 + - ^n ab + E ab = 0, (2.16) 

the (0, v)— field equation (see [7| for the explanation of this terminology), 

2 - h a h 9> b - h a b a hc , d h d c - n «* u c (cu d .j + 2 u d u f ) = q a , (2.17) 

the divergence of vorticity equation, 

u a , b h b a = u a u a , (2.18) 
and the magnetic part of the Weyl tensor, 

H ab = 2 u {a u b) - K h s b (co (t 9 ' c + a (t 9 ' c ) r] s)fgc u f . (2.19) 



Next the equations ( 2.11|) projected orthogonal to u a and along u a re- 



spectively give the equations of motion of matter, 

(jjl + p) u a + h ac (p )C + <K b C]b + q c ) + (u ab + a ab + ^6 h ab ) q b = 0, (2.20) 

and the energy conservation equation, 

fi + 6{fi + p) + 7c ab a ab + q a ;a + u a q a = . (2.21) 



Finally the various projections of the Bianchi identities fl2.12|) along and 
orthogonal to u a give analogous equations to Maxwell's equations for E ab and 
Haf, in ( 271 ). These consist of the div-E equation, 

h\ E° d .j h f d + 3uj s H h s - n ba ™ u a a d p H qd = ±h b c /i' c 

+ \ [-K M A + u b a cd n cd - 3 uj m q d + (T M q d -^9q b + n bd ii d } (2.22) 

the div-H equation, 

h» g H° d .j h f d -3u s E b s + n ba ™ u a a d p E qd = (ji+p) u b 

+\ V\ac W q a > c + \ r] b qac u* (u dc + a dc ) n a d , (2.23) 



the E -equation, 

h) E f9 ^ + h a b rf> rsd u r H a s . d - 2 Hf n t)drs u d u r 



-E? u b)s - 3 o- b)s + h tb E dp a dp + 9 E bt = --{pi + p) a 

-\ h th {fi + 9(fi + p)} - q( b ifi - i u«> q» - \ q ™ 
o 2 2 



tb 



+ -{u c{b + ct c(6 } u l) q c +-6u {t q h) ---h + tt c 



2 L 1 6 

2 O 



U ' u r 



(2.24) 



and the H -equation, 



h) H fg h] - h a b n t]Tsd u r E a s . d + 2 Ef rf )drs u d u r 

-Hf uo b)s - 3 Hf a b)s + h tb H dp a dp + 9 H bt = -q (t uo b) 

+ ^V {b radU t) U r {u Cd + a Cd }7T a c . 



(2.25) 



3 Perturbations of Isotropic Cosmologies 

We shall now assume that the space-time in section 2 is a perturbed 
RW space-time. Thus the background metric tensor g a t, is the Robertson- 
Walker metric, the background energy-momentum-stress tensor is Q2.2p spe- 
cialised to a perfect-fluid (by putting q a = = n ab ) with fluid 4-velocity 
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u a and the background Weyl tensor vanishes. The Ellis-Bruni [|lj approach 
to perturbations of this background is to work with gauge-invariant small 
quantities rather than small perturbations of the background metric. Such 
gauge-invariant quantities have the important property that they vanish in 
the background space-time. Thus for an isotropic background the Ellis- 
Bruni variables, which will henceforth be considered small of first order, 
are E a b, H a b, a ab , it a , u a b (or equivalently the vorticity vector u a ), X a = 
h b fx t b, Y a = h b a p b , Z a = h b 9 tb , n ab and q a . The equations satisfied by 
these quantities are given by fl2.15| )- (|2.20| ) and by (|2.22|) - (|2.25|) neglecting 



non-linear terms in these variables. Raychaudhuri's equation ( |2.14| ) and 



the energy conservation equation (|2.21|) are not immediately useful because 



they are not expressed in terms of these gauge-invariant variables. We can 
however get useful equations from (|2.14|) and ( |2.21| ) by calculating the spa- 



tial gradients of these equations and then retaining only linear terms in the 
gauge-invariant variables. This results in the two equations 



Z a + 0Z a -9u a -h ab (u a - a ) b + -X a + -Y a = , (3.1) 

' 2 2 

and 



X a + - 9 X a + (fx + p) Z a + 9Y a - fiii a + h ab (q c . c ) b = . (3.2) 
The background value of 9 to be substituted into ( |3.1| ) is given by 



= -±0*-IOi + 3p) , (3.3) 



which is got by specialising Raychaudhuri's equation ( |2.14 ) to the back- 
ground. In (|3T^ ) the background value of fi is given by 

fi = -9(ix + p), (3.4) 

which follows from Q2.21| ) specialised to the background. We shall assume an 
equation of state of the form p = p(fi) so that the Ellis-Bruni variables X a 
and Y a are related by 

Y a = ^X a . (3.5) 
an 

Finally we shall assume that in the background fx + p ^ so that we do 
indeed have a cosmological model in which the Einstein tensor determines a 
unique time-like 4-velocity of the matter. 

We look for solutions a ab , u a , uj a b (o- uj a ) , X a (and thus Y a by ( |3.5|) ), 
Z a , n ab and q a of the linearised versions of the equations ( |2.15| )-( |2.2(J| ) and 
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(|2.22 )- (|2.25|) along with fl3.1| ) and fl3.2| ) for which these variables depend 
upon an arbitrary function. This is because we expect that this dependence 
of perturbations will describe gravitational waves carrying arbitrary informa- 
tion. We note that E ab , H ab are derived variables. Specifically we assume 
that 

(Tab = s afe F(0) , ii b = a b F(0) , u ab = w ab F(<f)) oo a = w a F ((/>)) , 
X a = x a F((j)) , Z a = z a F((f)) , ir ab = Tl ab F((j)) , q a = Q a F{(f>) . (3.6) 

where F is an arbitrary real-valued function of its argument <fi(x a ). This form 
for the gauge-invariant variables was first introduced by Hogan and Ellis || 
where the perturbed matter distribution was taken to be a perfect fluid. As 
we mentioned in the introduction, the idea of introducing arbitrary functions 
into solutions of Einstein's equations describing gravitational waves goes back 
to pioneering work by Trautman M. We note that all of the quantities in 
( 3.6 ) are orthogonal to u a and that s a b , H a b are trace-free with respect to 
the background metric g ab (i.e. s a a = = n a a ). 

When ( |3.6j ) are substituted into the linearised versions of ( |2.15| )- fl2.20| ), 
(|2.22j )- (|2.25| ) and (j3~T|), ( pl2|) the following extensive but surveyable list of 
equations emerges as indicated. 

From the shear propagation equation: 

E ab = Q U ab + p ab ^j F + m ab F' , (3.7) 
with F' = dF/dcj) and 

11 2 

Pab = 0(o;b) + u (a a b) ~ ~8 U(a a b) — ~ O ;/ h ab — S ab — -9 S ab , (3.8) 
m ab = 0( a A 6 ) - - X Kb ~ 4> Sab ■ (3.9) 

Here and throughout A a = h\ 4> )b , x — 0,/ a ^ an d 4> — 4>,a u a . 
From the magnetic part of the Weyl tensor: 

H ab = q ab F + l ab F' , (3.10) 

with 

1 , 

q a b = W(a;b) +W( a U b ) - - 9W( a U b ) - S(/' C T] b ) fpc U J - h ab W C - c , (3.11) 
Lb = W [a \ b) - W c 0' C h ab - S (a P 11b)fpc U f (f)' C . (3.12) 
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From the div-E equation: 

m afe )fe = O, (3.13) 

n afe 0, b + p ab 0, fe + m afe ;6 = O , (3.14) 

W b , b +p ab , b = l -x a - l -eQ a . (3.15) 
From the div-H equation: 

r b 0, fe = O, (3.16) 

q ah <P, b + l ah ]b = \il a q hcU q Q b <jy c , (3.17) 

q ab -,b -^ + V)w a = \ 77V « 9 ^ • (3-18) 
From the E -equation: 

ri« + H u bt - - h bt Q a . a + - «(« g fe ) + - qm --9u^ g b ) = 

3 6 ' 2 ^ 2^ 6 

-P W - «r g ( 'W )rsd - 9p bt - l - (fi + p) s bt , (3.19) 

4>Yl bt - 1 - h bt Q a <P, a + l -Q {t \ b) = -4>p bt -m bt -6m bt -u r (q^ b s <P 4 +l {b s - 4 ) rf )rsd , 

(3.20) 

j)m bt + l {b s r 1 t > sd u r (i) 4 = Q . (3.21) 

From £/ie H -equation: 

q bt - u r p (b s , d rf> sd + 6q bt = 0, (3.22) 

g« + /« + g z w - Mr (p^ 0, d + m( fe s;d ) r^ sd = , (3.23) 

-mS h a rf> 8d u r <t>,d = . (3.24) 

From t/ie (0, v) -field equation: 

s ab < l> fi + ri ac *u c w d (l>j = 0, (3.25) 

g a = ^ a - A-^^c^;/ • (3.26) 

From the vorticity propagation equation: 

w a + ^ew a = l -r 1 acbd u c a b . 4 , (3.27) 



^w a = X -rf cM u c a h ^ A . (3.28) 



From the divergence of vorticity equation: 

<a = , 
W a ^, a = . 

From the spatial gradient of Raychaudhuri 's equation: 

z c + 9z c + ^x c + ^y c = h cb (a d ,^) & + 9 

<Pz c = a d . 4 \ c + h cb (a d <P, d ) b , 
(a d 0, d ) A c = . 

From the equations of motion of matter: 

y - + Q- = -^ + p ) a a -U ab , b , 
<PQ a = -U ab ^ b . 



c 



From the spatial gradient of the energy conservation equation: 

4 
3 



x c + - 9 x c + 9 y c + h cb (Q a . a ) b = -(fi + p)9a c -(fi + p) z c 



>x c + Q a ]a \ c + h b c (Q a >a ) 6 = O 
(Q a 0,a) A b = . 



On account of ( |3~1)| ) the variable y a appearing in ( |3.31| ), ( |3.34| ) and ( |3.36| ) is 
related to x a by 

dfx 

We must now examine the internal consistency of these equations. 
Our first enquiry parallels the work in ||. Putting 



ybt = m bt + i ft ^ 



we can write (|3.21|) and ( |3.24|) as a single complex equation 



2 4> V bt = % V trsd u r V h s 0,, + i V brsd u r V\ 0, d . 
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3.29) 
3.30) 



3.31) 

3.32) 
3.33) 

3.34) 
3.35) 

3.36) 

3.37) 
3.38) 



3.40) 



3.41' 



From this we calculate that 

^ bpql V bt u p = 2iV t [q \ l] . (3.42) 

When this is substituted into each term on the right hand side of ( |3.41|) we 
obtain 

2 V bt = 2 V bt + 2 4T V,d 0' d V bt . (3.43) 
Hence with <fi ^ and V ^ we must have 



. (3.44) 



The hypersurfaces 0(a; a ) = constant in the background isotropic cosmological 
model must be null (see section 4 below where the physical implications of 
(|3.44| ) are discussed). Thus A a = h b a 4> )b = )(1 + 4>u a ^ and so we now have, 
in addition to ( |3.44j ), the following vanishing scalar products: 

w a 0, a = O, a a <P A = 0, Q"0 ia = O, (3.45) 

on account of (|3.30|) , ( |3.33|) and ( |3.38|) . The latter two simplify ( |3.32| ) and 
(|3.37| ) respectively. 

We will next show that for consistency of our equations we must have 
Q a = 0. Let us, for convenience, write (|3.8|) as 



p ab = A ab a g ab ? (3^) 

3 

with 

A ab = a (a;b) + ^(a ^6) _ 1 q ^(a fl 6) _ ]_ Q f f h ab _ ^ ^ 

3 3 

Using the fact that in the background 9 >a = —0u a we have from fl3.46| ), 

p ab ]b = A ab b _-ab b _^ 9s ab b (g^g) 

In the background u a]b = | 6 h ab and so this equation can be written 

p ab b = A ab b _ ^b^ _ q g ab ^ _ (3 49) 

The Ricci identities satisfied by s ab give 

s ab ; Cb u c = (s ah . h ) , (3.50) 

and so ( |3.49| ) becomes 

p ab ;b = A ab , b -{s ab , b )-es ab , b . (3.51) 
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Now (|Q6| ) is 



„u6 



-Q a + - z a + A a 



with 



A a 



-r] acdf u c w d . f 



Putting flpD into (^5lD yields 



p ab . b = A ab . b + Q a + 9 Q a - -z a - -9 z a - A a - 9 A a 



Alternatively from Q3.15 ) and ( [3.34j ) we have 



1 



P ,b= -x a --9Q a + y a + Q a + (v + P) a a 
Thus ( |3.54| ) and ( |3.55| ) are consistent provided 



(3.52) 
(3.53) 

(3.54) 
(3.55) 



4 2 
-6Q a = - 
3 * 3 



z a + z a + - x c 



y a )+(li+p)a a + A a + 6A-A ab . ib . (3.56) 



Making use of (|3.28 ) to write w a in terms of a a and using the propagation 
equation ( |3.27| ) for w a along u a it follows that 



A a + 6A a 



1 1 

-d a + - 9 a° 

2 2 



1 

2 



--/i ac (a b . 



-(u-p)--9--9 2 \a° 
4 ^ ^ 6 9 f 



(3.57) 



Direct calculation from ( |3.47| ) yields 



A ab , b 



6 

Putting and Q335|) into 

equation (|3.3| ) gives 



-a a + ^a a + |i(/i-j9)-^--0 2 l a a 
2 2 l4 VP ^ 6 9 J 

+ -/i ac (a 6 . fe ) - -0 



^ . a ;fe w + - a . b . 



(3.58) 



and using the background Raychaudhuri 

\8Q a = \ [z a + 9z a +^x a + ^y a -h ac (a b , b ) c -9a a } . (3.59) 
It thus follows from (|3.31|) that 

9Q a = 0. (3.60) 

Since 9 > we must have 



Q a = . 



(3.61) 
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It now follows from ( |3.37| ) and (|3.38| ) that 

x a = , 

and so ( |3.39| ) yields 



(3.62) 
(3.63) 

(3.64) 



y a = o. 

Now (|3~32|) with ( ^45|) becomes 

4>z c = a d . d \ c , 

while ( |3.36| ) reduces to 

(n+p) (z c + 9a c )=0 . 
Contracting ( |3.64| ) with iC and using (3.44) results in 

2 c iC = <pa d ]d , 
and, noting (|3.45|) again, contracting ( p.65|) with iC gives 

(// + p) a d ;d = . 
With 0^0, /i + p^ we must have a d ; <2 = and so ( |3.64| ) becomes 

z c = . (3.68) 
Now (|3.65|) with 9>0andfi + p^0 yields 



(3.65) 



(3.66) 



(3.67) 



It now follows from Q3.28j ) that 



a c = . 



w a = . 



(3.69) 



(3.70) 



At this stage the only surviving gauge-invariant small quantities from the 
list ( |3.6|) are a ab and n ab or equivalently s ab , IP 6 . We see that now (|3.25|) 
and (|3.35|) require 

s ab 0,6 = 0, n ab 0,6 = 0. (3.71) 
Also (P^BD and (ET33) require s afc , II ab to satisfy 



s ab . b = , IT 







(3.72) 



We will write the remaining equations from ( |3.13|) - (|3.38|) , which have not 
reduced to = 0, in terms of s ab , IP 6 . We will then check that the equations 
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we obtain for s ab , IT ab (including (|3.71 ) and (|3.72|j ) are consistent. We begin 
with (F20|) with Q a = and substitute into it p ab from fl3j) with a a = 
and m ab from ( ^.9[ ) with a a = to obtain 



n 6t = 2 s w + - s bt + s 6t - u r (q( b s 6 >d + Z^ ;d ) rf rsd . 



With g a6 given by ( |3.11| ) with w a = we calculate 



-s"" <p tP u b + s tb ' p • 



+ 0S 



tb _ Jp;b , 



+ -96s 



tb 



(3.73) 



(3.74) 



Next using I in ( |3.12j ) with w a = we find 



r] brsd u r l 



st 



-0 S d t U b + S 6 t 0' d + S 6 t M d - S d t 0' 6 



(3.75) 



The first of ( |3.72| ) and the background expression u d - d = 9 help us to deduce 
from (|3~75|) that 



, !d s td u b - -6<j>s tb + s tb . d <f)' d + s tb < 
3 



(3.76) 



+0s* b + 0s' h + #0s* b 
Putting (|3.74j ) and ( |3.76|) together we get 

V brsd u r (q* a 6 4 + Z* g;d ) = 2 s*^ 6 4 + 0' d ;d S * b + 2 i* b + s' 6 + s' 6 , (3.77) 
which is symmetric in (£,&). Substituting this into (|3.73|) gives the equation 

1 ' (3.78) 



4 • iflo. s„, 



--0n tb , 



with s' tb := s t b ] d4>' d - This is a propagation equation for s t b along the null 
geodesies tangent to 6< d . 

We now turn our attention to (|3.23|) . Substituting for p ab , m ab from ( |3.8|) 
and (13.91) with a a = we see that 



p b s 6 d + m b s - d = -(6 d s b & 



9 <j>,d s b s -^9<j)u d s b s - s b s - d 



We notice in passing from this that now (|3.14|) , which reads 



p ao 6 >h + m ab . b = 



(3.79) 



(3.80) 



on account of the second of (|3.71| ), is automatically satisfied because of the 
first equation in ( p.71|) and in (|3.72| ). Writing out ( |3.23| ) with q ab substituted 
from (|3.11|) with w a = we have 



l bt + 9l bt - u r (0 s^ b S]d + p( b s 6 >d + m^ b s , d ) rf> sd = . 



(3.81) 
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Now using ( p.79| ) we can write this as 

l bt + 9l bt + u r {{<f>, d s^ b s y + 9 (f> >d s^ b s } ^ rsd = . (3.82) 

This can be re-arranged as 

(l bt + u r <j) 4 s (b s rf> sd )- + 6 (l bt + u r <j) 4 s {b s rf> sd ) = . (3.83) 

With l bt given in (|3.12| ) with w a = we see that ( |3.83| ) is identically satisfied. 

We next examine (|3.22|) . It is identically satisfied and this can be seen as 
follows: with p ab in (|3.8|) and a a = 0, we find 



P b s;d 



-s b s - d + ^eu d s b s -^0s b s;d 



(3.84) 



From the Ricci identities satisfied by s we have 

s b S ;d = {s b s -, d y +(-/< + -?) {u b s ds + s b d u s ) + -9 {s b s , d + s b s u d ) . (3.85) 
o 2 3 

This allows us to write 

V trSd UrP b s;d = -(V trSd U r S b s . id y-9 V trsd U r S b S]d . (3.86) 

When this is entered into (|3.22|) the equation can be re-arranged as 

(q bt + u r s^ h s . >d rf> sd y + 6 (q bt + u r s^ b S]d rf^ sd ) = , (3.87) 

which is an identity on account of (|3.11|) with w a = 0. 

We now consider ( J3.19Q with Q a = 0. If we first substitute for p ab into it 
from (|3.8|) with a a = we obtain 

■bt 



s bt - h -6s bi - 2 -es bt - 2 -6*s bt + l - fr+p) s bt + u r q ^ s;d rj 



(b t)rsd 



(3.i 



However using q in ( |3.11| ) with w a = we find that 

r, brsd u T q t s = -Us td u b + ^ 
and thus, using the first of ( |3.72| ) 



2 8 td u b + t_g s tb u d _ h td u b + h bt u d _ s td-,b + s tb-,d ^ ( 3 8Q ^ 



^2 tb + lg s tb + s tb e -tb 
9 3 



td;b I tb;d 



(3.90) 



The second last term here can be simplified using the Ricci identities satisfied 
by s ab and the first of ( |3.72| ) to read 



S td *- 4 



tb 



(3.91) 



15 



We see that now ( |3.90| ) is symmetric in (b, t) and on substitution into ( |3.88 
we arrive at a wave equation for s ab , namely, 



-9s ab 
3 



ab 



-tl ab --6U ab . (3.92) 



With w a = = Q a we have from (|318|) that 

q ab -,b = . 



(3.93) 



This equation is satisfied by q given by ( |3.11| ) with w a = 0. Substituting 



the latter expression for q ab into the left hand side of ( |3.93| ) and using (|3.91 
we have 



q 



1 

~4 

~ (If 1 ,,,. ■ 



p;cb & p;bc) V ^ ^/ 



gcb p 



+ R P9 c b s a9 ) v bcfp u f 



(3.94) 



using the Bianchi identities. Here R a gc b are the components of the Riemann 
tensor of the isotropic background. These components can be easily written 
in terms of the perfect-fluid energy-momentum-stress tensor of the back- 
ground because the background, being isotropic, is conformally flat. When 
this is done it is found that each of the Riemann tensor terms in (|3.94|) 
separately vanishes and so ( |3.93|) is satisfied. 

With the second of ( |3.72| ) holding we have from ( [3.15D 



P ab ;b 



. 



(3.95) 



That this is satisfied by p ab given by (3J3) with a a = is straightforward when 
one notes that s b s -b = {s b -b)' = 0, which follows from ( ^.85| ) after summation 
over (b, d) and the properties of s ab . 



Finally we must check ( p,17|) with Q a = 0. This reads 



q b<t>,d + l l 



b;d 



(3.96) 



To see that q ab and l ab satisfies this equation we start with ( p.77|) and multiply 
it by rjbpqt- Since the right hand side of ( |3.77| ) is symmetric in (t, b) we obtain 



( a ™ A) + p W - (a m 6 +l m ) 



q;mj ^p 



. 



(3.97) 



If this is multiplied by u p then ( |3.96| ) results. 

At this point all of the equations ( ^.13| )-( P.38| ) are satisfied provided 
s ab , U ab satisfy the algebraic relations with a in (|3.71|) , are divergence- 
free as indicated in ( |3.72| ) and satisfy the propagation equation fl3.78p for 
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s along the integral curves of <p ,a in the background space-time and satisfy 
the wave equation ( |3.92| ) . These equations reduce to those obtained in || 
when U ab = 0. We need to check that (|577Tp , (15772]) , ( gTTg ) and ( g/jg|) are 
consistent. 

Using the Bianchi identities we can show that 



(s ab ' d A b 



S ab ; b 



7 



■p 



n 



ab 



(u ab , b y + -eu ab lb 



(3.98) 
(3.99) 
(3.100) 



With the help of these equations it is straightforward to see that the wave 
equation (|3.92| ) is consistent with ( |3.72|) . Also using 



[s ab u b y d , d 



-0s ab . b , 
3 



one can easily see that the wave equation ( |3.92| ) is consistent with 



-,ab 



U ab u b . The wave equation (|3.92|) is also consistent with s ab cpj, 
This follows from 



u b 

o = n 



(3.101) 
= 



ab 



Q ab\d i 
S -d<P,b 



-2(s 



ab:d 



\d):b 



~ab 



;d),b , 



(3.102) 



which is obtained using the Bianchi identities satisfied by s ab and jCl . The 
propagation equation ( |3.78| ) allows us to write ( |3.102|) as 



S ;d<P,b 



■^0<P, b s ab + 



,U ab 



(3.103) 



One can also derive this by multiplying the wave equation ( |3.92|) by (j) tb and 
using the fact that (4>, b )' — 4>,b~ \ @ with, as always, = hl4>,c- 

The propagation equation ( |3.78| ) is clearly consistent with ( |3.71| ) because 
the integral curves of <f)' a are geodesies. That ( 3.78Q is also consistent with 

Ub-,c (p' c = \6 Ab. The consistency 



s ab u h 



= II u b follows from u' b 



of (|3.78|) with (|3.72|) requires the wave equation (|3.92|) . This is because on 
taking the divergence of ( |3.78| ) and using the Ricci identities one arrives at 



2 r 



W b , b + 



,n° 



> b s ab ' d - d + ^6s ab 



1 

3 



(3.104) 



Substituting for s ab ' d -d here from the wave equation (|3.92|) , this equation 
reduces to Tl ab ib = 0. 
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4 Pure Gravity Wave Perturbations 



As a result of the calculations outlined in section 3 the perturbations of 
the Weyl tensor now have "electric" and "magnetic" parts given by 

E ab = (lu ab + Pab )F + m ab F' , (4.1) 



with 



,2 

H ab = q ab F + l ab F', (4.2) 



2 

Pab = ~s ab --9s ab , m a b = -<f> s ab , (4.3) 

Qab = -S(a P ' ,C Vb)fpcU f , Lb = -S( a P Vb)fpcU f (P ,C ■ (4.4) 



Also s ab , II ab satisfy the consistent equations ( FT7ID , (gJg) , @7g ) and ( |3792|) 



and (p' a is a null vector field in the background RW space-time. From the 
first of ( ET7ID we see from Q and (4.4) that 



m a V, fe = 0, 1*^ = 0, (4.5) 



verifying that ( |3.13| ) and ( |3.16| ) are satisfied. Thus the F'-parts of F afe , H, 



ab 



above are type N in the Petrov classification with degenerate principal null 
direction ,ct . We therefore consider the F'— part of this perturbed field as 
describing gravitational waves having propagation direction ,ct in the RW 
background and the histories of the wave-fronts are the null hypersurfaces 
4>(x a ) = constant. This interpretation is based on the well-known analogy 
with electromagnetic radiation flO| . The F-parts of E ab , H ab are not in 



general type N and so do not necessarily describe gravitational waves. 

For the remainder of this paper we will consider pure type N perturba- 
tions (i.e. pure gravity wave perturbations) of the RW background. We could 
do this by requiring the F-parts of E ab , H ab in ( |4.1| ) and (4.2) to vanish. 
It is possible to exhibit solutions of our basic equations (3.71), (3.72), (3.78) 
and (3.92) having this property (see section 6 below) but it is less restrictive 
to require that the F-parts of E ab , H ab also be type N in the Petrov classi- 
fication with 0' a as degenerate principal null direction. This means that, in 
the light of the second of (|3.71|) , we should require 

p ab (j>, b = 0, <f b 0,6 = O, (4.6) 

with p ab , q ab given above in (|4.3j ) and (4.4). The first of these can be written 

s ab ct> Mc u c = Q , (4.7) 
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while the second gives us 



S ab <P,a' C -S aC <P, a ' b = . (4i 



To elucidate the meaning of ( |4.7|) and ( [4.8|) it is convenient to make use of 
a null tetrad in the background RW space-time. First we note that k a = 
—<fi~ 1 4 ) ,a and l a = u a — \k a are two real covariant vector fields satisfying 
k a k a = 0, l a l a = and k a l a = — 1. Let m a , m a be a complex covariant 
vector field and its complex conjugate (indicated by a bar) chosen so that 
they are null (m a m a = = rh a rh a ), are orthogonal to k a and l a and satisfy 
m a m a = 1. Now k a , l a , m a , m a constitute a null tetrad with respect to 
which s ab can be written (because s a a = , s ab Ub = , s ab kb = and so 
s ab l b = 0) 

s ab = sm a m b + sm a m b . (4.9) 

Thus \s\ 2 = \ s ab s ah . Substituting fl4~9| ) into (|4^ ) we easily see that ([L7|) is 
equivalent to 

s(p M m b l c = , (4.10) 
from which we conclude that provided s ^ (4^ s ab ^ 0) we must have 

M m b l c = O. (4.11) 

On using ( |4.9|) in (|4.8|) we find that, in addition to ( |4.11| ), 



s 4>,a-,b m a m = s 4>^b m a m . (4. 12) 

A simple way to satisfy this with s ^ is to require the null hypersurfaces 
<f)(x a ) = constant to satisfy 

<p^ b m a m b = . (4.13) 

This means that the complex shear of the null geodesic congruence tangent 
to <p' a in the background RW space-time vanishes. If we can find a family of 
null hypersurfaces cj){x a ) = constant in the background space-time satisfying 
( |4.11| ) and ([4.1 3|) then ( |4.6| ) will be satisfied. The solutions we then obtain 
of the equations (3.71), (3.72), (3.78) and (3.92) will be analogous to the 
Bateman waves H] of electromagnetic theory. 



5 Explicit Examples 

To exhibit explicit examples of pure type N perturbations of the RW 
space-times we first select in such space-times some naturally occuring shear- 
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free null hypersurfaces </>(x a ) = constant. We begin with the general Robertson- 
Walker line-element in standard form, 

^ = RW idZl)2 + (dx2)2 \ {dX ^-de, (5.1) 

O + i'- 2 ) 

where R(t) is the scale factor, r 2 = (x 1 ) 2 + (a; 2 ) 2 + (x 3 ) 2 and k — , ±1 is 
the Gaussian curvature of the space-like hypersurfaces t = constant. We can 
put these line-elements in the following interesting forms for our purposes 

rnj: 

ds 2 = R 2 (t){dx 2 +p 2 f 2 (dy 2 + dz 2 )}-dt 2 , (5.2) 

with po = 1 + ^ (y 2 + z 2 ) , K — constant , / = f(x). The following cases 
arise: (i) if k = +1 then K = +1 and f(x) = sinx; (ii) if k = then 
K = 0, +1 with /(x) = 1 when if = and /(x) = x when if = +1; (hi) if 
k = —1 then K = 0, ±1 with /(x) = | e 21 when K = 0, /(x) = sinhx when 
if = +1 and f(x) = coshx when K = —1. 

Case (i) above arises because when k — +1 the closed model universe with 
line-element ([5.1D has t = constant hypersurfaces with line-element which 
can be put in the form dl 2 = R 2 (t) ds^ with 

dsl = dx2 + sin2 x + sin2 ^ , (5-3) 

and we then use stereographic coordinates y, z such that y + iz = 2 e lip cot | 
in place of the polar angles ^ , v 9 - 

Case (ii) arises because when k = the open, spatially flat universe 
with line-element ( |5.1| ) has t = constant hypersurfaces with line-element 
dl 2 = R 2 (t) dsl where 

ds 2 = dx2 + d y 2 + dz2 > ( 5 - 4 ) 

or 

rfs 2 = rfx 2 + x 2 (dtf 2 + sin 2 d dip 2 ) , (5.5) 

and in the latter we introduce the stereographic coordinates y, z again in 
place of i? , </?. 

Case (iii) is due to the fact that in (|5.1| ) when k = — 1 the t = constant 
hypersurfaces can, modulo the factor R 2 (t), each be viewed as the future 
sheet of a unit time-like hypersphere Ti^ in four dimensional Minkowskian 
space-time M.^. Thus if the line-element of .M4 is written 

ds 2 Q = ( dzl ) 2 + (dz 2 ) 2 + (dz 3 ) 2 - (dz A f , (5.6) 

then Ti.3 is given by 

(z 1 ) 2 + (z 2 ) 2 + (z 3 ) 2 - (z A ) 2 = -1 , z 4 > . (5.7) 
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The different parts of case (iii) above are due to the different ways one can 
parametrise (|5.7 ). One possibility is with 

z 1 + iz 2 = (y + iz)pQ 1 sinhx , (5.8) 

3 - ( l (y 2 + z 2 ) - l) p^sinhx , (5.9) 



\4 

4 



z = coshx , (5.10) 



with p = 1 + My 2 + z 2 ). Substitution into (|5.6|) gives 



dsl = dx 2 + Pq 2 sinh 2 x (dy 2 + dz 2 ) . (5-11) 
The next possibility is 

z 1 +iz 2 = ^e x (y + iz) , (5.12) 

z 3 = ^e x [y 2 + z 2 - 1) + e~ x , (5.13) 

z i = ^e x (y 2 + z 2 + l)+e- x , (5.14) 



and now (R.6I) reads 



ds 2 = dx 2 + -e 2x {dy 2 + dz 2 ) . (5.15) 

Finally we can take 

z 1 + iz 2 = (y + iz) p^ 1 coshx , (5.16) 

z 3 = sinh a; , (5-17) 

z 4 = (^(y 2 + z 2 ) + l^j p^coshx , (5.18) 

with pq = 1 — | (y 2 + z 2 ). With this (|5T6|) takes the form 

dsl = dx 2 + Pq 2 cosh 2 x (dy 2 + dz 2 ) . (5.19) 



In the space-times with line-elements ( |5.2| ), with the special cases out- 
lined following fl5.2|), the hypersurfaces 

(f)(x a ) :=x- T(t) = constant , (5.20) 

with dT/dt = are null hypersurfaces. They are generated by null 
geodesies having expansion 

5^ = W + !' (5 ' 21) 
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Here f — df/dx , R = dR/dt. The integral curves of the vector field d/dt 
are the world-lines of the fluid particles. The components of this vector field 
are denoted by u a and using ( |5.2(J| ) we can show that 

2 4>,a;b = £a 4>,b + £b 4>,a + <j>J d 9ab , (5.22) 

with i 

ia = -J<P,a + R<P/u a . (5.23) 



With s ab (j) >b = = s ab u b we see on substituting ( |5.22|) that fl477|) and ([O 



are now satisfied. On account of ( |5.22j ) it follows that jQ is shear-free |12 



Alternatively we can easily verify (|4.11|) and (|4.13|) using the null tetrad 



described following ( [4.8| ) which is given via the 1 -forms 

k a dx a = Rdx — dt , l a dx a = — — (Rdx + dt) , 

m a dx a = -^Rpq 1 f (dy + idz) . (5.24) 
v2 

For convenience we have used the same coordinate labels {x, y, z, t} for 
all of the special cases included in (|5.2|). Of course the ranges of some of 
these coordinates will be different in the different cases (for example, in case 
(ii) x G (—00 , +00) if K — whereas x G [0, +00) if K = +1). Similarly 
the shear-free null hypersurfaces ( |5.20[ ) differ from case to case, and within 
cases (ii) and (iii), as can be seen by noting the intersections of these null 
hypersurfaces with the space-like hypersurfaces t = constant. In case (i) 
the intersection is a 2-sphere. In case (ii) it is a 2-sphere if K = +1 and a 
2-plane if K = 0. Thus ( |5.20 ) describes two quite different families of shear- 



free null hypersurfaces that can arise in an open, spatially flat universe. In 
case (iii) the intersection of ( |5.20|) with t = constant can be a 2-space of 
positive (K = +1), negative (K = —1) or zero (K = 0) curvature giving 
three different families of shear-free null hypersurfaces in a k = —1 open 
universe. A geometrical explanation for these subcases is given in |TT . 



We begin with <p(x a ) given by ( |5.20| ) and u a d/dx a = d/dt. Since s ab and 
IP 6 are orthogonal to u a and <p ,a and trace-free, with respect to the met- 
ric tensor given via the line-element (|5.2| ), each have only two independent 
components. If the coordinates are labelled x 1 = x , x 2 = y , x 3 = z , x 4 = t 
then the surviving components are s 33 = — s 22 = a(x,y,z,t) , s 23 = s 32 = 
f3(x,y,z,t) and IT 33 = -IT 22 = A(x,y,z,t) , IT 23 = II 32 = B(x,y,z,t). We 
can conveniently express these on the null tetrad ( |5.24| ). We have s ab given 
by (|3) with 

s = -R 2 p~ 2 f 2 (a + i f3) , (5.25) 
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and 
with 



n 



ab 



Iim a m b + Iim a m 



(5.26) 



n 



-R% 2 f(A + iB) . (5.27) 

Calculation of the first of ( p.72j ) to be satisfied by s ab shows that a , (3 must 
satisfy the Cauchy-Riemann equations 

^(Po 4 «)-^(PoV) = 0, (5.28) 

^(Po 4 /^) + ^(Po 4 «) = 0. (5.29) 
In addition IT afc and thus A , B must satisfy the same equations, 

^(Po 4 ^)-^(Po 4 5)=0, (5.30) 

^ {P " B) + lz {P " A) = - (5 ' 31) 
We shall find it convenient to work with a , (3q rather than a , (3 where 

a = afR 3 , (3 = pfR 3 . (5.32) 



Since / = f(x) , R = R(t) we have (|5.28| ) and (|5.29|) satisfied by ao and (3q 
and these equations can be written economically as 



Q 

^{Po 4 ("o + i@o)} = , 

with ( = y + i z, giving 

where Q is an analytic function of £. Now (|5.25|) reads 

s = -R-Y Q r 1 g(c,x,t) . 



(5.33) 



(5.34) 



(5.35) 



The propagation equation (13.78 ) for s ab along the integral curves of <p ,a gives 
A, B in terms of A)- Writing this in terms of II given by (|5.27 ) we find 
that 

n = -2 R- 2 plf-\DQ + RQ) . (5.36) 

Here the operator D is given by D = d/dx + Rd/dt = d/dx + d/dT with 
T(t) introduced in ( |5.20| ). Also the dot indicates differentiation of R(t). It 
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follows from this and ( [5.27 ) that A + i B is analytic in ( and so ( |5.30| ) and 
(|5.31| ) are automatically satisfied. The only remaining equation to satisfy is 
the wave equation ( |3.92| ). With s ab given by ( |4.9| ) and ( [5.35D and with IP 6 
given by (|5.26| ) and (|5.36| ) we find, after a lengthy calculation, that (|3.92|) 
reduces to the remarkably simple wave equation 

D 2 g + kg = o, (5.37) 

with k = 0, ±1 labelling the RW backgrounds with line-elements of the form 
(15.21). Thus we have for k = 0, 



g((, x, t) = a(C, x - T) (x + T) + 6(C, x - T) , (5.38) 

for k = +1, 

0(C, x, t) = a(C, x - T) sin (^—] + &(C, x — T) cos (^^) , (5.39) 



and for k — — 1, 

G(£,x,t) = a((,x-T) smh(^-^-) +b((,x-T) cosh f ^^'j , (5.40) 



where in each case a((,x — T), b((,x — T) are arbitrary functions. In deriving 
(|5.37| ) we have made use of the equations 

f" = -kf, (f) 2 + kf 2 = K, (5.41) 

which are satisfied in the cases (i)-(iii) described following (|5.2|) above. 

With s ab and U ab known we can calculate m ab , l ab , p ab and q ab in order 
to form the electric and magnetic parts of the perturbed Weyl tensor as 
indicated in ( |4.1|) and (4.2). We can write the result compactly as 

d 

E ab + iH ab = _ 2 R -2 p 2j-l ,g p j m a m b _ /g^x 

OX 

We emphasise that Q is given by ( |5.38| )-( |5.40| ) in the various cases and now 
F = F(x — T) so that F' = dF/dx. Also Po = 1 + f (y 2 + z 2 ) , / = f(x) 
described following (p.2|), and R(t) is the scale factor. It is immediately clear 
from ( 5.42 ) that the perturbations of the RW background which we have 
constructed here are pure gravity wave perturbations. We will discuss some 
of their properties in section 6. 
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6 Discussion 



The propagation equation (|3.78|) for s ab along the null geodesies tangent 
to <f) ,a shows that if s ab = then U ab = 0. An important converse property of 
the pure type N perturbations described in section 5 is that if U ab = then 
s ab = provided // + p^0. To see this we have from (|5.36 ) that IP 6 = is 
equivalent to 

Dg + kg = o. (6.i) 

Substituting this into the wave equation Q5.37D results in 

(R 2 - RR + k)G = . (6.2) 

For the background RW space-time the fluid proper density /i and isotropic 
pressure p satisfy 

^-{R 2 - RR + k) = fi + p , (6.3) 
R z 



as a consequence of Einstein's field equations. Hence we can rewrite (|6.2|) 
simply as 

(Lt + p)G = 0. (6.4) 

From this and ( |5.35|) it follows that s ab = provided /i + p^0. 

The perturbed Weyl tensor given via (|5.42|) for the pure type N pertur- 
bations can be infinite where po(y,z) is infinite (when y 2 + z 2 — ► +oo if 
K ^ 0) and where f(x) vanishes (when k = +1 at x — 0, when k — — 1 
at x — > — oo (K = 0) or at x = (K = 0)). There is one non-singular 
case corresponding to k = , K = +1 for which po = 1 and / = 1. In 
this case the expansion of the history of the wave-fronts ( |5.21| ) is entirely 
due to the expansion of the universe. This case is as close as one can get to 
plane-waves in the present context and is analogous to plane Bateman waves 
in electromagnetic theory. 

Had we wished to construct examples of type N perturbations for which 
the F-parts of (|4.1|) and (4.2) vanish we see from ( |5.42| ) that these would be 
given by ( |5.42|) with dQ/dx = 0. This condition would then be incorporated 
into the wave equation ( |5.37| ) and the appropriate solutions Q((,t) replacing 
(|5.38| )- (|5.40|) could easily be obtained. 



We note that we have used the assumption 6 > in the background 
cosmological models to conclude from (3.60) that Q a = 0. If the background 
were an Einstein static universe then 6 = and we would have Q a ^ 0. It 
is well-known (see, for example [fOp ) that the Einstein universe is unstable 
and it might be interesting to investigate how this instability manifests itself 
in the formalism we use in this paper. 
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There are exact cosmological solutions of Einstein's field equations known 



which contain gravitational waves (see [fL4}j , p5[ and references therein). 
These solutions describe universes with a stiff equation of state so that the 
speed of sound is equal to the speed of light. Our perturbations describ- 
ing gravitational waves propagating through isotropic cosmologies place no 
restriction on the equation of state of the isotropic background and thus 
we would not in general expect them to approximate to these known exact 
solutions. 
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